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We review recent lattice results for quark masses and low-energy hadronic parameters relevant 
for flavor physics. We do that by describing the FLAG initiative, with emphasis on its scope 
and rating criteria. The emerging picture is that while for light quantities a large number 
of computations using different approaches exist, and this increases the overall confidence on 
the final averages/estimates, in the heavy-light case the field is less advanced and, with the 
exception of decay constants, only a few computations are available. 

The precision reached for the light quantities is such that electromagnetic (EM) corrections, 
beyond the point-like approximation, are becoming relevant. We discuss recent computations 
of the spectrum based on direct simulations of QED-fQCD. We also present theoretical de¬ 
velopments for including EM effects in leptonic decays. 

We conclude describing recent results for the K —> tttt transition amplitudes and prospects 
for tackling hadronic decays on the lattice. 


1 Introduction and FLAG 

After its discovery in 2012, the Higgs boson was believed to provide a portal to New Physics. 
This is even somehow assumed when formulating the hierarchy problem of the Standard Model 
(SM). However, this far, all measurements of the Higgs boson properties lie within 20% of the 
SM expectations, as reported by ATLAS ^ and CMS^. Instead, there is a number of 2-3 sigmas 
tensions in rare processes (see for example^), with the most prominent examples being in the 
angular analysis of the —)• K*^pt^pL~ decay ^ and in the enhancement of the B —)• D^*^Th'r 
decays^ . Significances depend on treatment of several non-perturbative effects. Extrapolating to 
the future, (some of) these rare decays won’t be so rare anymore. Belle 2 will report results from 
about 2018 and coexist with the LHC and High Luminosity (HL)-LHC, after Long Shutdown 
3 in 2023-2025. Progress on the theoretical side is needed in many instances, to match the 
expected experimental accuracy. 

There are many different groups, all over the world, using different lattice methods, that 
calculate hadronic matrix elements relevant for a number of weak decay processes of K, and 
mesons. With so many groups calculating similar matrix elements using different methods, 
and all providing phenomenologically relevant results with complete error budgets, it is useful to 



try to produce global averages/estimates and to review virtues and shortcomings of the different 
computations in a transparent way, which should be accessible also to the non-experts. This is 
the goal of the FLAG initiative. 

1.1 The FLAG review 

The Flavor Lattice Averaging Group started its activity in 2010 focusing on light-quark quan¬ 
tities and providing averages from lattice results with comprehensive error budgets®. A second 
similar initiative was started at around the same time^, focusing on both heavy- and light-quark 
quantities. The two groups joined for the second edition of the FLAG-review® (FLAG-2). One 
of the main goals of FLAG is to assess the reliability of systematic error estimates, in particular 
concerning continuum extrapolations, chiral extrapolations, finite volume effects and renormal¬ 
ization. This is done through quality criteria by assigning to each computation a symbol for 
each one of the systematics above. For example, the symbols and the criteria adopted for the 
light-quark quantities in FLAG-2 are: 

• Chiral extrapolation: 

★ < 200 MeV 

o 200 ’ MeV < < 400 MeV 

■ 400 MeV < 

in addition it is assumed that the chiral extrapolation is done using at least three points. 

• Continuum extrapolation: 

★ 3 or more lattice spacings, at least 2 points below 0.1 fm 

o 2 or more lattice spacings, at least 1 point below 0.1 fm 

■ otherwise 

in addition it is assumed that the action is 0(a)-improved (i.e. the discretization errors 

vanish quadratically with the lattice spacing). 

• Finite-volume effects: 

★ > 4 or at least 3 volumes 

o > 3 and at least 2 volumes 

■ otherwise. 

• Renormalization (where applicable): 

★ non-perturbative 

o 1-loop perturbation theory or higher with a reasonable estimate of truncation errors 

■ otherwise. 

For heavy-light quantities the criteria are similar, with some additional ones concerning dis¬ 
cretization effects and treatment of heavy quarks. In general criteria are expected to change in 
time and possibly become stricter as lattice computations reach new levels of accuracy. In the 
end, all the published (in journals) results with no red symbols enter the final estimates/averages. 
In some cases, the averaging procedure leads to results which in the opinion of the authors do 
not cover all uncertainties. In these cases, in order to stay on the conservative side, averages are 
replaced by estimates (or ranges), which are considered fair assessments of the current knowledge 
acquired on the lattice. These estimates are based on a critical (and to some extent subjective) 
analysis of the available information. 

In detail, the FLAG-2 collaboration counted 28 members representing the major lattice 
groups in the world. Different Working Groups were in charge of reviewing different sets of 
quantities: Quark masses (WGl), Vus,Vud (WG2), yPT Low Energy Constants (WG3), Bk 
(WG4), foyy Bb (WG5, our group), R(s)) D semileptonic and radiative decays (WG6), 

and finally as (WG6). In the following we will focus on the subset of quantities presented during 
the talk. A more recent update of lattice results concerning heavy-light quantities can be found 
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Figure 1 - Lattice results for the strange quark mass ms and the average up and down quark mass in the MS 
scheme at the 2 GeV scale. The bottom panels represent non-lattice (and PDG) results. The FLAG-2 final 
estimates, including filled green points only, are given by the grey bands. Figure from®. 


in®. In Fig. 1 we show the results for the strange quark mass and the average up and down quark 
mass. That also serves the purpose of clarifying the difference between averages and estimates 
discussed above. Indeed, in the Nf = 2 + 1 case an error has been included in the final estimates 
accounting for the quenching of the charm quark (see® for details). 

A second instructive example is taken from the Vus, Vud working group. These CKM matrix 
element can be extracted from leptonic as well as semileptonic decays and therefore the WG2 
focuses on kaon and pion decay constants as well as on form-factors relevant for the K —)• 
Tiiv transition. In particular the form-factor /-|-(0) at zero momentum transfer is relevant for 
phenomenology and for comparisons to xPT. In Fig. 2 we show the summary plots for these 
quantities from ®. The results can be used to check the first row unitarity of the CKM matrix in 
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Figure 2 - Lattice results for /+(0) and the ratio of decay constant fn/fk- The FLAG-2 final estimates, including 
filled green points only, are given by the grey bands. The blue points are non-lattice estimates. Figure from®. 

the SM. Neglecting Vub, the Nf = 2 -|- 1 estimates give |I4dp -|- |I4sp = 0.987(10). In addition, 
as discussed in the review, the consistency of leptonic and semi-leptonic determinations of |I4s| 
is a check of the equality of the Fermi constant describing interactions among leptons, and the 
one describing interactions among leptons and quarks. This gives an important constraint on 
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possible modifications and extensions of the SM. 

As mentioned, the current situation in the heavy-light sector is much less satisfactory. While 
some quantities like decay constants have been computed by a number of collaborations using 
a large variety of methods, for more complicated ones like form-factors, even for the “simplest” 
pseudoscalar to pseudoscalar, tree-level induced, semileptonic transitions, only a few determi¬ 
nations exist. In Fig. 3 we show the FLAG-2 summaries for /b(„) and for the B —)• Triu form 
factor /_|_(g^). For the latter only two computations, based on the same ensembles of config¬ 
urations (but using different treatments of heavy quarks) exist. The situation is similar for 
other quantities like the Bb^^^ mixing parameters (see®). Very much like experimental results, 
the confidence increases when several results from different collaborations/experiments become 
available. Significant progress in this directions is indeed expected from the lattice community 
in the next few years and will be visible already in the next FLAG review (expected for 2016). 
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Figure 3 - Lattice results for (left) and for the B —>■ nlv form factor /+(q^) (right). The BCL 

parameterization® is used for the latter. The form-factor is expressed as a function of the z variable, obtained 
from g® through a conformal transformation depending on a real parameter topt- Figure from*. 


2 Inclusion of EM interactions in lattice QCD computations 

Most of the lattice calculations concerning the properties of the light mesons are performed in 
the isospin limit of QGD and neglecting EM interactions. However, at the precision reached 
(e.g., the FLAG-2 estimates for the pion and kaon decay constants have an error < 1%), they 
cannot be ignored anymore. For example, the EM corrections to the mass of the charged pions 
are estimated to be 4 - 5 MeV. The current approach mostly relies on yPT for correcting lattice 
data in order to include both EM and strong isospin breaking effects. Obviously it would be 
desirable to deal with the corresponding terms directly at the level of the simulations. 

The BMW Collaboration reported in about the first direct simulations of QED and QCD 
with four non-degenerate flavors, in a fully dynamical formulation. The goal is the first-principle 
computation of the neutron-proton mass difference, a tiny (0.14%) effect, which is crucial in 
explaining the Universe as we know it. This impressive computation is summarized in Fig. 4, 
where the results for the contour lines of the neutron-proton mass difference are given in terms 
of the niu — md = Sm splitting and the EM coupling a (both normalized to their physical 
value) “. Within the same approach the authors of also compute the mass splittings in the 
E, E, Ecc and D channels. They also provide an estimate of the numerical cost for such a 

“The separation among EM and strong (QCD) isospin breaking effects is ambiguous by 0{a5m). The pre¬ 
scription adopted in the LO (in isospin breaking) calculation in^'^ fixes the EM correction to the mass difference 
m^- — m-^+ to vanish. 












































computation. Considering the various extrapolations/interpolations in masses and couplings, 
the poor statistical signal for small values of the electromagnetic constant and the need for very 
large volumes, such a calculation turned out to be 300 times more expensive then their pure 
QCD computation of the spectrum of stable hadrons in the theory with two dynamical flavors. 



Figure 4 - Contour lines for the neutron-proton mass difference resnlting from a direct QED-fQCD computation 
with l-l-l-l-l-l-l (i.e., non-degenerate) dynamical flavors. Figure from^^. 

Indeed finite volume effects are one of the main issue in simulating QED on a lattice because 
of the long-range nature of the EM interactions. In particular in a hnite volume with periodic 
boundary conditions (in space) zero modes of the gauge field exist, which can not be eliminated 
through (standard) gauge-fixing conditions. In^^ the finite-volume zero mode is removed through 
a non-local constraint. In fact, a rigorous, all-order, proof of the renormalizability of the theory 
in this setup is still missing. An alternative could be to give a mass to the photon perhaps a 
la Stueckelberg (and references therein). The massless limit (which would have to be taken 
numerically) is smooth in this case, at least in the continuum. On a similar line, it may be 
interesting to reconsider soft covariant gauges, as proposed and studied in for non-Abelian 
gauge theories. 

Additional issues, due to infrared divergences, must be taken into account when trying to in¬ 
clude QED corrections into the computation of transition amplitudes. These are already present 
in the case of the decay constants (or better, the case of leptonic decays), as discussed in^^. Let 
us consider the widths describing the 7r+ —>■ decay at 0(q:) and label them as Tj, with i the 
number of photons in the final state. It is well known that to obtain physical quantities radiative 
corrections from virtual and real photons must be combined. Therefore, at this order we are 
interested in Tq and ri(AE), where the energy of the photon in the final state, and in the rest 
frame of the vr'*', is integrated from 0 to AE. For the sake of the argument on which the approach 
in is based, it is sufficient to look at the subset of diagrams shown in Fig. 5, all contributing 
to Tq. The first (from the left) gives the pure QCD contribution and it is factorizable into an 
hadronic part (encoded in the matrix element of the axial current between the vacuum and a , 
i.e., the decay constant) and a leptonic one, because a W boson is exchanged in between the two 
vertices. The second one is again factorizable and could be viewed as an ©(a) correction to the 
decay constant, however it is infrared divergent. These divergences are removed by considering 
diagrams as the rightmost one and diagrams where a photon is emitted either from a quark 
line or the lepton line (diagrams contributing to Ti). But the third diagram, where a photon 
is exchanged between a quark and the lepton, is not factorizable, so there is really not much 
physical sense in “QED corrections to decay constants”, rather one should consider corrections 
to the whole transition process. 

In principle, both Tq and ri(AE) could be computed on the lattice, however the latter would 
be computationally very expensive. Instead, the authors of propose to use the pointlike {pt in 
the following formulae) approximation to calculate ri(A£'). Values of AE around 10 - 20 MeV 
are experimentally accessible and for such soft photons the coupling to hadrons is conceivably 





Figure 5 - Pure QCD (left) and examples of factorizable (middle) and non-factorizable (right) contributions to 
Fo at 0(a), see text. Figure from^"^. 


well described by the pointlike approximatior/’. In order to ensure an accurate cancellation of 
the infrared divergences, and since Tq is computed on the lattice through numerical simula¬ 
tions, whereas Ti{AE) is computed in perturbation theory (and in the pt-approximation), it is 
convenient to introduce an intermediate step and re-write 

r{AE) = ro + ri(A.B) = {ro-r(;*} + {rP* + rf(A£:)} 

= ^i^m {ro(L)-rg‘(L)} + {< + rf(Aii;)} , (i) 

where L is the linear extent of the lattice. As pointed out in^^, the small momenta contributions 
to ro(L) and rg*(L) are the same, hence the infrared divergences cancel in the difference. The 
same is true for the infinite volume combination Tg* -|- T^^{AE), therefore the two terms in the 
brackets on the r.h.s of eq. 1 are separately infrared finite (and, incidentally, gauge invariant) 
and have a well defined infinite volume limit. 

The correlation functions needed in the lattice computation of ro(T) and in particular the 
three-point functions required for the non-factorizable terms are explicitly constructed in 
The implementation of the method is computationally demanding, but seems within reach of 
present resources, and first studies should soon be performed. 

3 Results and perspectives for hadronic decays on the lattice 

Many phenomenologically interesting transitions involve hadronic two-body final states and 
the lattice would be extremely useful in providing hrst-principle computations which would 
serve in clarifying existing tensions (e.g., the one mentioned in the angular analysis of —)• 
K*^{K it) by LHCb), or give an ab-initio explanations of long-standing puzzles such as the 
AJ = 1/2 rule and the value of e'/e in K ^ tttt decays. However, there is no simple relation 
among Euclidean correlators and the desired Minkowski-space transition matrix elements, a fact 
which is known as the “Maiani-Testa no-go theorem” A solution to this problem, for the 
case where one two-particle state only (say, tttt) is kinematically accessible or coupled to the 
initial state {K), was developed by Liischer and Lellouch in a series of papers I6,i7,i8,i9^ ^ 

first step a relation is established, in Minkowski-space, between the finite volume dependence of 
the energy levels of two-particle states (vrTr) and the infinite-volume S-matrix and phase shifts. 
Since energy levels are directly computable in Euclidean-space, this allows to measure elements 
of the S-matrix on the lattice. The kaon is introduced in a second step and it is coupled to the 
two-pion states through a small, perturbative, Weak-Hamiltonian term Riy- The lattice volume 

*’The Lagrangian describing the interaction of a pointlike meson with the leptons is non-renormalizable by 
power counting, very much like the chiral Lagrangian, which could indeed have also been used here (although 
more complicated from the analytical point of view when considered in a finite volnme). For the process described, 
however, these interactions are inserted at tree-level only and therefore there is no need for additional connterterms. 
The only requirement for the method to work is that the contributions from small momenta are the same in the 
full theory and in its approximation. Still, for larger values of AE some approximations may be more accurate 
than others. 




has to be tuned such that one of the two-pion energy levels gets degenerate with the kaon, that is 
what is usually called “matching the kinematics”. At this point degenerate perturbation theory 
can be used, and as in the first step, a relation (in terms of “Euclidean” quantities) is established 
among the perturbative corrections to the two-particle energy levels in finite volume and the 
perturbative corrections to the infinite-volume S-matrix, which is to say a relation among the 
finite- and infinite-volume versions of the {'KTi\H\[r\K) matrix element. The latter then gives the 
K —7- Trvr transition amplitude. 

The approach has been generalized to the case of multiple strongly-coupled decay 

channels into two scalar particles and to the case of external currents injecting arbitrary four- 
momentum as well as angular momentum. These are first steps towards lattice computations of 
amplitudes for processes such as D —)■ Trvr and D ^ KK and towards study of meson decays as 
BO K*^{K7r)n+n-. 

In the case of the K —)■ tttt transitions, numerical results became recently available with good 
control over all the systematics including continuum limit extrapolations. These are outstanding 
results of many years of efforts and attempts. In^^ the amplitude A 2 for a kaon to decay into 
two pions with isospin 1 = 2 has been computed on two lattices with resolutions a = 0.11 fm and 
a = 0.084 fm respectively. The calculations have been performed using 2+1 flavors of domain 
wall fermions with pions at the physical mass and L ~ 5 fm. The matrix elements of three 
different operators have to be combined in this case and the final result, extrapolated to the 
continuum limit, reads 

ReAs = 1.50(4)stat(14)syst x 10“® GeV , (2) 

ImAs = -6.99(20)stat(84)syst x 10-^® GeV , (3) 

which is well consistent with both the very accurate (but different) experimental numbers for 
ReA2 from charged and neutral kaon decays. The error on the lattice value is dominated 
by systematics, in particular by the uncertainty in the perturbative evaluation of the Wilson 
coefficients, currently known at NLO. 

The computation of the Aq amplitude is much more demanding as 10 operators including 
QGD penguins producing quark disconnected diagrams, need to be considered'^. However, after 
the conference, two independent preliminary results appeared, both using a single lattice 
spacing and both reporting on a computation of the A/ = 1/2 amplitude Aq. 

4 Conclusions 

Flavor Physics is still playing a prominent role in the indirect search for New Physics. At the 
same time, and while finalizing the analysis of LHC run I data, new signals from direct searches 
may emerge (as for example in the search for resonances presented in^®, and interpreted within 
composite dynamics models in^®), which will hopefully be confirmed by run H. 

The picture provided here is obviously incomplete and the result of our taste and interests, 
but we hope to have given a flavor of the important role, the main challenges and the exciting 
future directions and perspectives for lattice gauge theories within Flavor Physics. As the 
keywords seem to be precise and rare, the lattice community is tackling all subleading effects 
(e.g., isospin breaking) and theoretical obstructions (e.g. multi-hadron decay channels) to give 
an indispensable contribution to the quest for New Physics. 
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